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Abstract 

This work is devoted in the derivation of novel upper and lower bounds for the Rice /e-function. 

These bounds are expressed in closed-form and are shown to be quite tight. This is particularly 
evident by the fact that for a certain range of parameter values, the derived lower bound virtually 
behaves as a remarkably accurate approximation. As a result, the offered expressions can be considered 
useful mathematical tools that can be efficiently employed in various analytical studies related to 
natural sciences and engineering. To this effect, they can be sufficiently applied in the area of digital 
communications over fading channels for the derivation of explicit representations for vital performance 
measures such as bit and symbol error probability, among others. 

Index Terms 

Closed-form representations. Rice /e-function, lower and upper performance bounds, special functions 

I. Introduction 

Over the last decades, special functions have played a crucial role in studies related to natural 
sciences and engineering. In the wide field of wireless communications, their application in 
various analytical studies often allows the derivation of closed-form expressions for critical 
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performance measures such as probability of error, channel capacity and higher order statistics. 
Furthermore, the majority of special functions are found to be built-in in widely used mathemat¬ 
ical software packages such as Maple, Mathematica and Matlab. As a result, the algebraic 
representation of any associated relationships as well as their computational realization, have 
been simplified enormously. 

Among other, a special function related to the area of telecommunications is the Rice Ie- 
function. This is denoted as Ie(k,x ) and was firstly proposed by S. O. Rice about six decades 
ago ffl. I n the field of digital communications it has been largely utilized in the study of zero 
crossing rates. In addition, it has extensively appeared in the analysis of angle modulation systems 
as well as in radar pulse detection and error rates in differentially encoded systems, If2l-lf5l. 

The definition of the Rice /e-function is typically given in an integral form with the involved 
integral being finite and comprising an exponential function and a modified Bessel function of 
the first kind and order zero. Alternative algebraic representations include two series expressions 
which were reported in 0. These series are infinite and are expressed in terms of the modified 
Struve function and the modified Bessel function of the first kind, f6), 0. Moreover, the author 
in [5j| reported useful closed-form identities that relate the Je-function with the Marcum Q- 
function, Il8l- lfl6l and the references therein. 

However, to the best of the authors’ knowledge the Rice Je-function is not explicitly expressed 
in terms of other special and/or elementary functions. In addition, it is not included as a built-in 
function in the aforementioned popular software packages. 

Motivated by these statements, this work is devoted in deriving upper and lower bounds 
for the Rice Je-function. These bounds are expressed in closed-form and are shown to be 
particularly tight. This is evident by the fact that for a certain range of parametric values, 
the lower bound virtually behaves as a remarkably accurate approximation. Furthermore, the 
algebraic form of the derived bounds is relatively simple and therefore, convenient to handle 
both analytically and numerically. This feature renders them useful for meaningful utilization 
in analytical studies relating to the area of digital communications over fading environments 
E2-EQ and the references therein. 

The remainder of this paper is organised as follows: Section II revisits the definition of the 
I e (k , x) function along with its existing alternative representations. Novel upper and lower bounds 
are derived in Section III while their performance is analysed in Section IV. Finally, closing 
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remarks are given in Section V. 


II. The Rice /c-function 
The Rice /e-function is defined as a, a, 

pX 

Ie(k,x) = / e~ t Io(kt)dt (1) 

Jo 

which is valid for 0 < k < 1. The notation J 0 (.) stands for the modified Bessel function of the 
first kind and zero order, 0, 0. An equivalent integral representation to (Q]) was proposed in 
0, namely, 


Ie(k, x ) = 


g— x{l—cosQ) 


dd 


( 2 ) 


\/l - k 2 n Jo 1 - kcosd 

In the same context, alternative representations for the Ie(k,x ) function were reported in 0- 
These representations are expressed in infinite series form and are given by, 


Ie[k , x) = 


xi r 


X 


X 


n T,2n 


2y/l ~ k 2 ^ n\ 2 n \fl - F Wl - k 2 n+5 


-- L n+ i (xy/1 - fc 2 ) + L n _ 1 (W1 - A; 2 ) 


(3) 


and 


Ie(k,x ) =xe x -^ 


f x (l-k 2 ) 


n =0 




n+1 


f n+1 




+ xe 


/o(/CX ) + 2 


n=0 


r(" + |) 

a; (1 — k 2 ) 


‘n+1 


(/ex) 


7 r (n + I) 


(4) 


where the notations L n (.), /„(.) and T(.) denote, the modified Struve function, the modified 
Bessel function of the first kind and the Euler’s gamma function, respectively 0, 0. According 
to 0, the Je-function can be computed by using the above two series alternatively. To this effect, 
equation 0 converges relatively quickly for the case that the term x\Jl — k 2 is large while the 
term kx is small. On the contrary, equation @ converges quickly when the term x\/l — k 2 is 
small while the term kx is large. 

Nevertheless, it is noted that this method for computing the Rice /e-function is quite inefficient. 
The reasons underlying this statement are the following: i) two relationships, instead of one, 
























4 


are required for the computation; ii) the above series are relatively unstable due to their infinite 
form; in) the L n (.) function is not a built-in function in popular mathematical software packages 
such as Maple , Mathematica and Matlab. 

An adequate way to resolve was reported in |J5]|. There, the Rice Je-function is related to the 
Marcum Q-function by the following two expressions 

Ie(k, x ) = . 1 [2 Q(a, b ) - e~ x I 0 {kx ) - l] (5) 

yl — At 

and 


Ie(k , x) = 1 [Q(a, b) - Q(b, a)] 

yl — k z 

where, 


( 6 ) 


a = 



(7) 


and 


b = \fx\j 1 — Vl — k 2 (8) 

The notation Q(a, b) denotes the Marcum Q-function of the first order which is defined as, 0, 

EQl 


/ cxj 2 2 

xe I 0 (ax)dx (9) 

III. Bounds to the Rice Je-FUNCTiON 

It is recalled here that the authors in [11] established the monotonicity criteria for the Q m (a, b) 
function and subsequently derived tight upper and lower bounds to it. By following a similar 
procedure in this Section, novel upper and lower bounds for the Rice Je-function are derived in 
closed-form. To this end, it is noted that a primary requirement is the derivation of an alternative 
representation of Je-function. 
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A. An alternative representation of the Rice Ie-function 

Lemma 1 . For x > 0 and 0 < k < 1, the following relationship holds, 


rx 

Ie(k,x ) = 1 — e~ x Io(kx ) + k e^IifkFjdt (10) 

Jo 

Proof. By integrating once by parts ([]]). it follows immediately that 


Ie(k, x) = 


e dt 


L«/o 


Io(kt) 


>o Uo 


e dt 


Io(kt) 


dt 


(ID 


According to the basic principles of integration, the first integral in (fTTI) is evaluated straight¬ 
forwardly as, 


e- f dt = -e~ x (12) 

Regarding the derivative of the modified Bessel function of the first kind and order n, it is 
recalled that according to flU, 0, it follows that 

4-I n (kx) = ^ [I n _i(kx) + I n+ i(kx)] (13) 

dx 2 

which for x = t and n = 0 yields, 

ffl°(kt) = kli(kt) (14) 

Subsequently, by substituting (fl2l) and (fl4l) in flTTI) and recalling that by definition /o(0) = 1, 
equation (fTOl) is finally deduced and thus, the proof is completed. ■ 

B. An upper bound for the Rice Ie-function 



Theorem 1 . For x > 0 and 0 < k < 1, the following inequality holds, 


Ie(k,x ) < 1 — e X l 0 (kx ) + 



er/(cy / x) 

c 



( 15 ) 


where c = \Jl — k and d = f l + k. 
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Fig. 1. The monotonically decreasing behaviour of I n (x) w.r.t n. 

Proof. It is recalled that the modified Bessel function of the first kind is strictly decreasing with 
respect to its order n. This is clearly illustrated in Figure 1. To this effect, it is easily shown that 
the inequality I n+a (x ) < J„(x) holds for a > 0. Based on this, the following inequality applies 
to (flQl) . 



( 16 ) 


Notably, for the special case that n + 0.5 6 N, the I n (x) function can be equivalently expressed 
in closed-form according to []6l eq. (8.467)], namely, 



(17) 


Therefore, for the special case where n — 0, it follows that 



( 18 ) 
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Evidently, the derivation of a closed-form upper bound to Ie(k,x ) is subject to evaluation in 
closed-form of the integral in (fl6l) . To this end, with the aid of (fl8l) it follows that 


/ e t Ii(kt)dt = 

Jo 2 

which after long algebraic manipulations yields 


/V 

i— 

»se 

1 

1 

1 

L e 

v/27 rkt 


dt 


(19) 


ex It, 

-i- ■' ’■ 


e Ii(kt)dt = y- 
where c = \/l — k and d = \/l + k. Also, 


erf(cy/x) erf{dy/x) 
c d 


( 20 ) 


erf(x ) = -^= [ e t2 dt 

^ Jo 


( 21 ) 


denotes the error function 0. Finally, by substituting (l20l) into (fl6l) . equation ( 031 ) is 



Fig. 2. The monotonically increasing behaviour of Q m (a, 6) w.r.t m. 


deduced, which completes the proof. ■ 


Remark. The authors in IfTO derived closed-form bounds for the generalised Marcum Q-function, 


















Q m (a, b ). Therefore, at a first sight it appears that by making the necessary change of variables 
and substituting accordingly in © and/or I©, an alternative expression for (fl5l) can be deduced. 
Nevertheless, it is noted here that the corresponding resulted expression is clearly less compact 
and elegant compared to equation (IT5l) . 

C. A lower bound for the Rice I e-function 


Theorem 2. For x > 0 and 0 < k < 1, the following inequality holds, 


Ie(k, x ) > 


2 Q{b + a) + 2 Q(b — a) — e x Io(kx ) — 1 

yT 


( 22 ) 


where 



Fig. 3. Behaviour of the bounds in eqs. (13) and (19) for k = 0.5 


i roo 2 

Q(x ) = / e~~dt 

nJ x 

denotes the one dimensional Gaussian Q-function [©. 


( 23 ) 
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Proof. According to the aforementioned monotonicity property of the I n (x) function, it im¬ 
mediately follows that that I six) < Ifx). As a result, by making the necessary change of 
variables and substituting in equation (fTOl) . the following closed-form inequality is deduced, 



Fig. 4. Behaviour of the bounds in eqs. (13) and (19) for x = 7 


fX 

Ie(k,x ) > 1 — e~ x I 0 (kx ) +k e~ t l3(kt)dt (24) 

Jo 2 

A similar inequality can be also deduced with the aid of equations © and ©. To this effect, 
it is firstly recalled that the Marcum Q-function is strictly increasing with respect to the order 
m, which is clearly illustrated in Figure 2. Based on this, the following inequality is obtained 
straightforwardly, 


Qi(a , b ) > Qi (a, b ) 

Subsequently, by substituting (1251) into ©, it follows that 


(25) 


Ie(k, x) > 


1 


VT^ 


2Qi(a, b) — e X l 0 (kx ) — 1 


(26) 
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It is recalled that based on the monotonicity criteria of the Marcum Q-function, the authors 
in [11] derived corresponding upper and lower bounds for Q m (a,b) along with a closed-form 
expression for the special case that m is an odd multiple of 1/2, i.e. m + 0.5 £ N. In this context, 
they thoroughly showed that 


Q i (a, b ) = Q(b + a) + Q(b - a) (27) 

To this effect, by substituting (1271) into (l26l) yields (l22l) . which completes the proof. ■ 

Remark. A lower bound for the Ie(k, x ) function could be theoretically derived by following the 
same methodology as in Theorem 1. Nevertheless, this approach ultimately leads to an analytic 
expression that is divergent. 



Fig. 5. Behaviour of the lower bound in (19) for x = 40 


IV. Numerical Results 

This section investigates the tightness and overall behaviour of the derived closed-form bounds. 
To this end, Figure 3 depicts the bounds in (fl5l) and (l22l) with respect to x, for k = 0.5 and 
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Fig. 6. Absolute relative error, e ar , between Ie(k,x) and (19) for x = 80. 

in comparison with results obtained from numerical integrations. Evidently, it is shown that the 
upper bound is tighter than the lower bound for small values of x. However, it is noted that 
as x increases, its tightness degrades while the lower bound becomes significantly tighter. As 
illustrated in Figure 4, this behaviour also holds for relatively small values of x as well as for 
different values of k. 

Interestingly enough, for large values of x-typically x > 40- the lower bound in (l22l) nearly 
becomes an exact representation for Ie{k,x) as it appears to behave as a remarkably accurate 
approximation. This is illustrated in Figure 5 where for the case that x = 80, equation (l22l) 
provides an excellent match to the corresponding theoretical results. This is also justified by the 
low level of the involved absolute relative error which is given by e ar —\ Ie(k,x ) — eg A l22l) | 
/Ie(k,x), and is depicted in Figure 6. 

V. Closing Remarks 

This work was devoted in the derivation of novel upper and lower bounds for the Rice Ie- 
function. These bounds are expressed in closed-form and are shown to be quite tight. This is 
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particularly the case for the lower bound since for large values of the argument x, it becomes 
a remarkably accurate approximation to /e-function over the whole range of values of k. As a 
result, it virtually behaves as an exact expression. Another advantage of the derived bounds is that 
their algebraic representation is relatively simple and therefore rather tractable both analytically 
and numerically. This is undoubtedly quite useful since the Rice Je-function is not included as 
built-in function in widely mathematical software packages such as Maple, Mathematica and 
Matlab. As a result, the offered expressions can be considered useful mathematical tools that 
can be meaningfully utilized in analytical studies related to the area of digital communications 
over fading channels. To this end, they can be indicatively applied efficiently in the derivation 
of explicit expressions for critical performance measures such as the bit and symbol error 
probability, among others. 
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